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Abstract 


Background: The molecular basis of severe acute respiratory syndrome (SARS) coronavirus 
(CoV) induced pathology is still largely unclear. Many SARS patients suffer respiratory distress 
brought on by interstitial infiltration and frequently show peripheral blood lymphopenia and 
occasional leucopenia. One possible cause of this could be interstitial inflammation, following a 
localized host response. In this study, we therefore examine the immune response of SARS-CoV 
in human peripheral blood mononuclear cells (PBMCs) over the first 24 hours. 


Methods: PBMCs from normal healthy donors were inoculated in vitro with SARS-CoV and the 
viral replication kinetics was studied by real-time quantitative assays. SARS-CoV specific gene 
expression changes were examined by high-density oligonucleotide array analysis. 


Results: We observed that SARS-CoV was capable of infecting and replicating in PBMCs and the 
kinetics of viral replication was variable among the donors. SARS-CoV antibody binding assays 
indicated that SARS specific antibodies inhibited SARS-CoV viral replication. Array data showed 
monocyte-macrophage cell activation, coagulation pathway upregulation and cytokine production 
together with lung trafficking chemokines such as IL8 and ILI 7, possibly activated through the TLR9 
signaling pathway; that mimicked clinical features of the disease. 


Conclusions: The identification of human blood mononuclear cells as a direct target of SARS-CoV 
in the model system described here provides a new insight into disease pathology and a tool for 
investigating the host response and mechanisms of pathogenesis. 


Background novel coronavirus [1-3] with genome sequence revealing 
The causative agent for SARS has been identified as a _no strong homology to existing known coronaviruses [4- 
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6]. Coronaviruses belong to the family of enveloped 
viruses called Coronaviridae, and have the largest known 
single-stranded viral RNA genomes (27 to 32 kb). Coro- 
naviruses, have both "early" and "late" phases of gene 
expression. Regulatory proteins are synthesized as "early" 
non-structural proteins, while the structural proteins are 
synthesized as "late" proteins. "Late" structural proteins 
are usually required in greater amounts thus, there is a 
necessity to regulate the expression of the viral genes 
quantitatively. After the viral entry via endocytosis or 
through specific receptors, the 5'-end of the viral genome 
is translated directly giving rise to twenty-three viral pro- 
teins, including the RNA dependent RNA polymerase 
(RdRp), and other functional products involved in tran- 
scription, replication, viral assembly and cell death. Coro- 
naviruses can be classified into species and three major 
antigenic groups based on, serology, natural hosts, mono- 
clonal antibody recognition and nucleotide sequencing 
[7]. Most coronaviruses have restricted host ranges as they 
infect only one host species or, at most, a few related spe- 
cies, they are an important group of animal pathogens. 
Group one (I) includes human coronavirus 229E (HCoV), 
porcine transmissible gastro-enteritis virus (TGEV) and 
feline enteric coronavirus (FECoV). Group two (II) 
includes bovine coronavirus (BCoV), murine hepatitis 
virus (MHV), and HCoV-OC43; and Group three (III) 
includes avian infectious bronchitis virus (IBV) [7]. Some 
coronaviruses like HCoV have restricted tissue tropism, 
including macrophages [8], although most strains that 
infect humans cause only mild respiratory infections. 


However, SARS has rapidly caused a world-wide problem. 
The earliest known cases of SARS was reported in Guan- 
dong Province, China in November 2002, becoming 
more widespread by March 2003, when it was introduced 
to Canada, Singapore, Taiwan and Vietnam via Hong 
Kong. The largest number of infected patients has been in 
China with a worldwide incidence totalling more than 
8,400 by July 2003. Infection by the virus induces high 
morbidity and mortality, the latter being estimated at 
15% by the World Health Organisation. SARS is character- 
ized by high fever, non-productive cough or dyspnea and 
in many cases may progress to generalized, interstitial 
infiltrates in the lung, thus needing intubation and 
mechanical ventilation [2]. The characteristic compres- 
sion of alveolar sacs seen in atypical pneumonia is largely 
due to fluid build up outside the alveoli. One possible 
cause of this could be interstitial inflammation, following 
a localised host response. To date, the details of the host 
response to SARS-CoV infection is still largely unknown 
and consequently the most appropriate treatment regime 
remains to be established. Typically a pro-inflammatory 
cytokine profile (up regulated TNFa, I11, IL6 and IFNs) is 
seen in viral infections such as influenza [9], together with 
perhaps limited amounts of IL8 and other chemokines 
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[10] that may depend on which cell type is infected [11]. 
In experimental systems the immediate innate immune 
response has been shown to be directed by the monocyte- 
macrophage-dendritic lineage to a range of different 
organisms [12,13] and consists of a core set of pathways 
common to all, together with pathogen specific pathways. 
This data points to critical time points in the response, 
with the first 12 hours representing primary events while 
longer periods the consequence of this activity and a sec- 
ondary (perhaps larger) cascade of responses. 


We postulated that the pulmonary damage in SARS may 
not be a direct effect of the virus on the alveoli, but repre- 
sents a secondary effect of cytokines or other factors prox- 
imal to but not from the lung tissue mediated by the host 
either as the primary or secondary response [2,14]. In this 
study, we have addressed this question by developing a 
human in vitro model system that will in the future allow 
detailed investigations of the host response to be made. 


Methods 

Cell culture and virus infection 

PBMCs were obtained by Ficoll-Hypaque separation of 
whole blood. 2 x 105 PBMCs were seeded into each well 
of a 24-well culture plate, 0.5 ml of complete RPMI-1640 
(Life Technologies-Invitrogen, USA) added to each and 
cultured overnight at 37°C (5% CO,). A seed stock of 
SARS-CoV (strain SIN 2774) passaged in Vero E6 cells was 
used for infection. Vero E6 culture supernatants were 
added to each well in 50 pl volume at a concentration of 
0.1 or 0.01 MOI (based on plaque forming units) and a 
control plate (media only). Each culture was set up in 
duplicate. After 4 hours of incubation, one set of the 
duplicate wells of the control plate and a 4 hours incuba- 
tion plate were harvested while the rest received 0.5 ml 
media top-up and incubated for a further 2, 4, 6 and 8 
days. 


Cell harvest and RNA isolation 

Harvesting was performed by gently flushing the wells 
with a Pasteur pipette to removed non-adherent cells, fol- 
lowed by a rinse of 1 ml RPMI. The rinsed fraction was 
pooled with the first harvest aliquot and spun at 1500 
rpm. The cell pellet was washed twice with 2 ml RPMI to 
remove Virus in the supernatant. 1.5 m1 Trizol (Invitrogen, 
USA) was then added to the adherent cell fraction as well 
as the non-adherent cell fraction to lyse cells and stabilize 
the RNA. Extraction of total RNA was then performed fol- 
lowing manufacturer's protocol and the resultant RNA 
dissolved in 40 ul water. 


Real-time quantitative polymerase chain reaction (PCR) 

The amount of SARS-CoV in each cell fraction was meas- 
ured by real-time quantitative PCR assay. 2 nl of RNA was 
reverse transcribed and amplified in 20 ul using 0.9 uM 
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each of forward (5'-GGTTGGGAT TATCCAAAATGTGA-3') 
and reverse (5'-AGAACAAGAGAGGCCATTATCCTAAG- 
3') primers, and 0.25 uM of TaqMan® MGB probe: 5'-(6- 
FAM)AGAGCCATGCCTAACAT(NFQ)-3') in a one step 
PCR using master mix from Applied Biosystems (USA) 
according to manufacturers' recommendations. Reactions 
were performed using an ABI PRISM 7900 sequence detec- 
tion system (48°C for 30 min, followed by 95°C for 10 
min and 40 cycles of 95°C for 15 sec and 60°C for 1 min) 
and quantitation achieved using standard curves gener- 
ated from in vitro transcribed RNA. 


High density oligonucleotide array hybridization 

At each time point (4 hours, 8 hours, 12 hours, 24 hours), 
5 x 10? cells of mock-infected and infected cells were har- 
vested and lysed using Trizol (Invitrogen, USA). Total 
RNA was isolated according to the manufacturer's recom- 
mendation. Quality of the total RNA was judged from the 
ratio between 28S and 18S RNA after agarose gel electro- 
phoresis. 20 ug of total RNA was labeled with Cy-3 or Cy- 
5 using the Superscript II reverse-transcription kit (Invitro- 
gen, USA) and hybridization was carried out overnight 
(16 hours) at 42°C on high-density oligonucleotide 
arrays (~19,200 gene features, Compugen) using univer- 
sal human reference (Stratagene, USA) as a reference. 
Hybridized arrays were scanned at 5 um resolution on a 
GenePix 4000A scanner (Axon Instruments) with variable 
photo-multiplier tube voltage to obtain maximal signal 
intensities, and the resulting images were analyzed via 
GenePix Pro v4.0 (Axon Instruments) as described in the 
manual. 


Microarray data analysis 

Raw data were analyzed on GenePix analysis software ver- 
sion 4.0 (Axon Instruments) and uploaded to a relational 
database. The logarithmic expression ratio for a spot on 
each array was normalized by subtracting the median log- 
arithmic ratio for the same array. Data were filtered to 
exclude spots with a size of less than 25 um and any poor 
quality or missing spots. Since the correlation of the over- 
all data from reciprocal labeling was good, values 
obtained from reciprocal labeling experiments were aver- 
aged. In addition, the data were distilled to the set of gene 
features that were present at all 4 time points in both the 
viral infected samples and the negative controls. The 
results were represented as the logarithmic ratio of gene 
expression between the viral infected samples and their 
corresponding negative controls at the various time 
points. Application of these filters resulted in the inclu- 
sion of ~12,900 of the total ~19,200 gene features in sub- 
sequent analyses. 


To discover patterns of gene expression, the values associ- 
ated with each gene feature f were translated so that their 
means were zero. Similar genes, whose translated gene- 
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features exhibited same induction-repression pattern, 
were grouped together. Genes g;, | were said to be similar 
if they satisfied the following condition: 


Y Usir * gi.) = N, where U(x) = 1 if x > 0, U (x) = 0 


otherwise, 


where g;, and g;, denote the translated values of gene fea- 
tures g;, g; at time t respectively; and N is the number of 
time points for which the expression of a gene was 
observed. Similar genes, based on the above criteria, are 
grouped together. Within each group, the genes were 
ordered in the descending order of their expression range 
(defined as the difference between the maximum and 
minimum ratios of gene expression). This algorithm is a 
special case of the Friendly Neighbor algorithm currently 
under development. The final plots were generated using 
the original expression ratios while preserving the cluster- 
ing and ordering discovered by the above algorithm. To 
determine whether a gene observed to be responsive 
could appear merely by chance, 100,000 expression pro- 
files were generated by randomly sampling the expression 
ratios from the entire dataset with replacement. The P 
value of a gene is the fraction of the random profiles 
whose logarithmic expression range is as good as, or better 
than that of the selected gene. 


Results and Discussion 

Kinetics of SARS-CoV replication 

We obtained PBMCs from 6 healthy volunteers by Ficoll- 
Hypaque separation of whole blood. Of the 6 donor 
PBMCs tested, all were able to support SARS-CoV replica- 
tion when infected with multiplicity of infection (MOI) of 
0.1. The first sampling, taken from cells infected for 4 
hours, showed an average copy number of 32 x 103 (Fig. 
1A) and represents the initial inoculum level. Over the 
course of the next 8 days, there was a steady rise in viral 
load, reaching as high as 480 x 103 copies per well in one 
donor, which could only be explained by active replica- 
tion of the SARS-CoV intracellularly. This work is sup- 
ported by recent in vivo evidence suggesting that SARS- 
CoV may have infected and replicated within PBMCs of 
SARS patients [15] and cells from humans and animals 
[16]. An indication of the PBMCs lineage involvement 
was provided by repeating the experiment using the 
monocyte-macrophage cell line THP-1 [17], in which viral 
replication was similar to the primary cell culture over the 
first 4 days (Fig. 1B). In the primary cultures, the non- 
adherent cell fraction which comprises mainly lym- 
phocytes and granulocytes showed dramatically less viral 
replication in our assay as did all cells infected at MOI of 
0.01 (Fig. 1A). 


The kinetics of viral replication was variable among the 6 
donors (Fig. 1C). There was a lag phase of 2 days in the 
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Figure | 


Kinetics of SARS-CoV replication. A. Viral copy number per well during culture. Two MOI infecting doses were used, 0.1 
(4,) and 0.0! (C,A) with 2 x 105 PBMCs/well. Harvested PBMCs were separated into two fractions, adherent cells (mostly 
macrophages), (I, A) and suspended cells (mostly lymphocytes), (A,L]). Points represent the mean of six subjects, each done 
in duplicate. Increasing viral copy number is clearly seen in adherent cells using an MOI of 0.1. There was little increase in viral 
load in non-adherent cells infected with SARS- CoV. B. Replication of SARS-CoV in the THP-| cell line. Using an MOI of 0.1, 
fold increase from day 0 in SARS-CoV infected THP-I cells (#) was similar to the PBMCs adherent fraction (A). Infections 
were done in duplicate and assayed in duplicate. Error bars represent SEM. C. Variation in the PBMC response to SARS-CoV 
infection. Using an MOI of 0.1 adherent cells (Ml) but not suspended cells (A) showed differing lag phase length and viral repli- 
cation in 6 subjects (a to f). Points show means of two replicates. 
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case of donors a, c and e; and 4 days for donors b, d and f 
before any significant increase could be detected. The viral 
replication generally peaked at either day 4 or day 6. The 
exception was donor b, in which the virus seemed to rep- 
licate at a much slower pace compared to the other 5 
donor samples. Equally interesting was the different levels 
of virus attained. Donor d seems to stand out from the 
rest, reaching a peak of 480 x 103 copies per well which is 
4 times more than that attained by donor e, with 120 x 
103 copies per well. Such variation strongly suggests that 
there is an underlying host-pathogen interaction influenc- 
ing the kinetics of SARS-CoV replication efficiency. These 
in vitro observations may reflect the wide range of patient 
outcomes after SARS-CoV infection [18]. 


Antibody blocking experiments were also performed in 
which SARS-CoV was pre-incubated with convalescent 
patient sera for 30 minutes before introduction to the 
PBMCs and after a 4 day incubation period, the adherent 
cell fraction was harvested and assayed for SARS-CoV viral 
titer. Results clearly showed that even at high dilution, 
convalescent sera inhibited SARS viral replication (data 
not shown), presumably by blocking viral entry. This sup- 
ports other reports indicating that SARS-CoV is not endo- 
cytosed through antibody mediated mechanisms and 
confirms a protective role for antibodies elicited either by 
the infection or through immunization [19,20]. 


SARS-CoV specific gene expression changes 

To further elucidate the molecular processes of SARS-CoV 
infection, PBMCs from 3 healthy individuals were 
infected separately in vitro with SARS-CoV (0.1 MOI) and 
harvested at 4 hours, 8 hours, 12 hours and 24 hours time 
intervals post-infection. As controls, uninfected aliquots 
of the same PBMCs were also harvested at the correspond- 
ing time points. Total RNA extracted from the PBMCs of 
the 3 individuals were pooled, labeled and hybridized to 
human oligonucleotide arrays consisting of ~19,200 gene 
features. Reciprocal dye swap replicate hybridizations 
were performed to minimize technical noise. Analysis of 
variance in expression levels for each gene across all the 
time points indicated the ~1200 genes which showed the 
largest variability (Fig. 2A and 2B). 


In order to focus the analysis, we queried the entire data 
set for genes related to the immune response by keyword 
searches on their gene ontology descriptions with the aim 
of describing the specific host-pathogen interaction. In 
common with other studies of respiratory pathogens [9- 
13], our data points towards two critical time points in the 
response, with the first 12 hours representing a primary 
pro-inflammatory cytokine profile while longer periods 
represent the consequence of this activity and a secondary 
cascade of responses [9-13]. We observed that within the 
first 12 hours of SARS-CoV infection, evidence of this 
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monocyte-macrophage activation was seen, indicated by 
enhanced expression of CD14, TLR9 plus NFKB1 and 
GATA signaling (Fig. 2C and Table 1). In addition, the 
MRC2 endocytotic receptor was upregulated as was the 
complement pathway (Clq, C6). Taken together, these 
data suggest an early activation of the innate immunity 
pathway. This activation was accompanied by an unusual 
cytokine transcriptional profile (Fig. 2C and Table 1). 
While IL1B (up regulated for the first 12 hours) would be 
expected following macrophage activation [21], TNFa, 
IFNy and IL6 were noted by their surprisingly low level of 
expression. This is in spite of the presence of elevated IL19 
which is thought to enhance their up regulation [22]. In 
some clinical investigation, concentrations of TNF and IL 
6 measured during active disease were found to be rela- 
tively low [23,24], reflecting our findings. This paper did 
not report on IFN levels, however, we found them to be 
low (Supplementary figure [see Additional file 1]). This is 
of particular interest as IFNs have been shown to have sig- 
nificant anti-SARS-CoV effects [25]. Such effects suggest 
that alteration of the IFN response and perhaps other 
immune modulators might provide opportunity for novel 
treatment and management regimes for SARS patients to 
be developed. 


A number of CC chemokines (CCL4, CCL20, CCL22, 
CCL25, CCL27) and their receptors (CCR4 and CCR7) 
were highly expressed in response to the infection (Fig. 2C 
and Table 1), indicating a rapid mobilization and 
increased trafficking, in particular of the monocyte-mac- 
rophage lineage very early on in the infection [26]. CXC 
chemokines (CXCL9, CXCL12) were also highly expressed 
suggesting significant increase in neutrophil homing as 
well. These are likely to be lung directed as IL8 and IL17 
were also highly expressed [27-32]. Specific trafficking of 
these cells to the lung may account for the localized nature 
of the response [33]. 


Surprisingly, a number of blood coagulation genes were 
highly expressed early during our in vitro infection (Fig. 
2C and Table 1), in particular TBXAS, which has been 
implicated in vasoconstriction, platelet aggregation, 
membrane lysis and increased permeability [34,35]; 
fibrin (FGB and FGG) and the coagulation pathway 
directly (SERPINs D1 and A3 together with Factors 10, 3 
and 2). This gives a pro-coagulation profile, which mimics 
the clinical-pathological observations: at autopsy, many 
SARS patients have unusually disseminated small vessel 
thromboses in the lungs without evidence of dissemi- 
nated intravascular coagulation [1,36]. Again, these 
expression profiles provide an experimental framework to 
explore an important aspect of SARS pathobiology and 
treatment. 


Page 5 of 11 


(page number not for citation purposes) 


BMC Infectious Diseases 2004, 4:34 http://www.biomedcentral.com/1471-2334/4/34 


Variance 


0 2000 4000 6000 8000 10000 12000 14000 


No. of gene features 


Time(h) _ Representative example of 
Cc Clusterno. Cluster profiles 4 8 12 24 — genesineach cluster 


1 


B Time(h) 


BNIP Hs.77572 


CXCL9 Hs.77367 


CCR4 Hs.184926 


FGB Hs.7645 


IF144 Hs.82316 


CXCL12 Hs.237356 


TLRO Hs.87968 


TBXAS1 Hs.2001 


IL19 Hs.71979 


IL2Ry Hs.84 


SERPINI2 Hs.158308 


IL4 Hs.73917 


CXCL3 Hs.89690 


SERPIND1 Hs.1478 


14 


SreerssaCheiit 
TTT LL 


Time (h) 04.812 


y 
x 
w 
3S 


| 


Figure 2 

SARS specific gene expression changes. A. Plot showing the variance of expression ratio for each of the gene features across all 
the time points. Dotted line indicates the ~1200 most variable gene features. B. Expression ratios of ~|200 gene features as 
grouped by the algorithm as discussed in Methods. Rows represent individual genes, columns represent individual time points. 
Each cell in the matrix represents the mean expression level from 3 subjects for a gene feature at a particular time point (non- 
infected PBMCs responses have been subtracted from infected responses). The red and green color bars reflect high and low 
expression levels respectively, while black indicates equivalent expression level. The magnitude of the log-transformed ratio is 
reflected by the degree of color saturation. C. Levels of PBMCs mRNA expression of 1087 immune-related genes at 4 hours, 8 
hours, 12 hours, and 24 hours in response to SARS-CoV infection were grouped and ordered using the algorithm described in 
Methods. Rows represent individual genes, columns represent individual time points. Each cell in the matrix represents the 
mean expression level from 3 subjects for a gene feature at a particular time point (non-infected PBMCs responses have been 
subtracted from infected responses). The red and green color bars reflect high and low expression levels respectively, while 
black indicates equivalent expression level. The magnitude of the log-transformed ratio is reflected by the degree of color sat- 
uration. The line graph indicates the average expression ratios for each group. The area above the axis indicates upregulation, 
while the area under the axis means downregulation. 


Page 6 of 11 


(page number not for citation purposes) 


BMC Infectious Diseases 2004, 4:34 


Table |: Representation of selected immune-related genes upregulated during the first 12 hours post-infection. 
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Gene Unique Gene ID. 


Function 


Expression range (log,) 


Cluster no. 


Immune Response Genes: 


TGFBRI Hs.220 
CD36 Hs.75613 
TLRO Hs.87968 
MST Hs.349110 
IL 19 Hs.71979 
CIQR Hs.97199 
C6 Hs. 1282 
GATA3 Hs.169946 
MMD Hs.79889 
CCL 25 Hs.50404 
CXCL 12 Hs.237356 
CCL 27 Hs.225948 
C4BPB Hs.99886 
CCR 4 Hs.184926 
CCL4 Hs.75703 
IL 17C Hs.1 10040 
IL 1p Hs.126256 
MRC 2 Hs.7835 
CXCL9 Hs.77367 
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Table |: Representation of selected immune-related genes upregulated during the first 12 hours post-infection. (Continued) 


CUL 4 Hs.183874 Tumour suppressor 5.35 6 
BAF 53A Hs.274350 Cell-cell signaling 5.26 6 
MAPK | Hs.324473 Apoptosis induction 5.11 14 
CUL 2 Hs.82919 Tumour suppressor 5.09 6 
P2RX | Hs.41735 Apoptosis, ion channel activity 4.43 6 
TYROBP Hs.9963 Intracellular receptor signaling 2.51 13 
CNIH Hs.201673 Intracellular signaling 2.13 9 
ADRalA Hs.52931 Cell-cell signaling 2.08 3 
DAPK | Hs.153924 Apoptosis induction 2.01 9 
Defense-related Genes: 

TRF | Hs.350470 Defense response, maintenance 7.65 9 
TFF 3 Hs.82961 Extracellular defense response 6.98 7 
GAGE | Hs.12823 1 Cellular defense response 6.82 6 
DEFB3 Hs.283082 Extracellular defense response 6.01 5 
RU 2 Hs.61345 Extracellular defense response 5.83 6 
DEFa4 Hs.2582 Extracellular defense response 5.81 6 
Interferons-related Genes: 

IFNaRI Hs.1513 JAK-STAT cascade, recptor 5.41 8 
IRF 7 Hs.166120 Transcription regulation 5.05 6 
IF 144 Hs.82316 Invasive growth response 2.69 5 
IRF 6 Hs.1 1801 Transcription regulation 2.06 14 
IL 18 Hs.83077 Angiogenesis 1.76 6 
ISGF3y Hs.1706 Regulatory 1.28 9 
IRF 5 Hs.334450 Transcription regulation 1.23 9 
IFNaR2 Hs.86958 Receptor activity 1.06 2 
IFNyRI Hs. 180866 Receptor activity 0.80 5 
IFNy Hs.856 Growth regulation 0.53 4 
IRF 3 Hs.75254 Transcription regulation 0.40 5 
IRF 2 Hs.83795 Transcription regulation 0.32 4 
IRF | Hs.80645 Transcription regulation 0.30 3 


Order is arranged from the highest to lowest value based on the expression range (log,). 


It is interesting to note that the TLR9 was highly expressed 
in comparison to other TLR receptors, implying some 
degree of TLR specificity for the virus (Fig. 3A). TLR9 is 
known to respond to CpG signaling motifs (GICGTT) 
[37-39] and one possibility is that the virus is activating 
directly through this mechanism. In support of this, we 
found that the SARS-CoV viral sequence contains the 
highest number (7 copies) of such specific signaling 
motifs compared to other coronaviruses and significantly 
more than several other viruses involved in respiratory 
diseases (Fig. 3B). It is conceivable that TLR9 may be aid- 
ing host recognition of the virus via the CpG groups and 
contributing to the initiation of the innate host inflamma- 
tory response. An alternative explanation is that TLR9 is 
being stimulated by mechanisms unrelated to CpG 
recognition. 


The emerging picture from this study implicates a central 
role for the immune response in the pathobiology of a 
SARS infection. While detailed in vivo studies of the host 
response are now required, the in vitro model described 
here will allow responses to specific modulators (such as 
therapeutics) to be investigated. In future developments 


of the model, it will be interesting to compare the host 
response to different SARS-CoV isolates with inactivated 
preparations of the virus. In other diseases, in vitro models 
have revealed a number of key processes relevant to the 
clinical diseases [9,12,13] and it is likely that the 
responses identified here will prove to be equally impor- 
tant. Although some clinical parameters have now been 
used as prognostic markers [40-42], further study of the 
regulatory mechanisms for chemokine-cytokine 
production will likely improve their accuracy and perhaps 
allow development of new treatment protocols. 
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Figure 3 

Expression of TLR9 in response to SARS-CoV infection. A. Expression range (log,) for TLR9Y, TLR2 and TLR4. The expression 
range for TLR9 was greater than expected (* represents a P-value for TLR9Y of 0.016). B. Comparison of the CpG motif 
(GTCGTT) copy number in coronaviruses and other viruses linked to respiratory diseases. Accession numbers are as follows: 
SARS coronavirus SIN2500 — AY283794, Human coronavirus 229E — NC_002645, Murine hepatitis virus - NC_001846, Avian 
infectious bronchitis virus — NC_001451, Bovine coronavirus — NC_003045, Human rhinovirus B — NC_001490, Human 
parainfluenza virus | - NC_003461, Human respiratory syncytial virus - NC_001781, Human metapneumovirus — 
NC_004148. 
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Additional material 


Additional File 1 


List of all immune related genes after SARS-CoV infection Compre- 
hensive list of 1087 immune related genes that were altered in PBMCs in 
response to SARS-CoV infection at 4 hours, 8 hours, 12 hours, and 24 
hours. Genes were grouped and ordered using the algorithm described in 
Methods. Rows represent individual genes, columns represent individual 
time points. Each cell in the matrix represents the mean expression level 
from 3 subjects for a gene feature at a particular time point (non-infected 
PBMCs responses have been subtracted from infected responses). The red 
and green color bars reflect high and low expression levels respectively, 
while black indicates equivalent expression level. The magnitude of the 
log-transformed ratio is reflected by the degree of color saturation. The line 
graph indicates the average expression ratios for each group. The area 
above the axis indicates upregulation, while the area under the axis means 
downregulation. 

Click here for file 
[http://www.biomedcentral.com/content/supplementary/1471- 
2334-4-34-S1.pdf] 
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